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1.  Historical  Introduction 


The  simplest  way  of  defining  a  caterpillar  tree,  Pn  (mi,  m2,  ...,  mn),  is  through 
the  concept  of  the  derivative  of  the  graph. ^  Thus  when  all  the  end  points  of  a  graph 
C  are  deleted  another  graph  G'  results  called  the  derivative  of  G.  A  caterpillar  tree 
is  defined  to  be  a  tree  graph  the  derivative  of  which  is  a  path.  The  name  caterpillars 
was  suggested  by  A.  Hobbs. ^  Thus  a  caterpillar  tree  Pn(m-|,  m2,  ...,  mn)  may  be 
constructed  by  the  addition  of  m-|  monovalent  vertices  to  the  the  first  vertex,  vi  of 
path,  Pn,  m2  monovalent  vertices  to  V2  of  Pn  and  so  on.  An  example  of  a  caterpillar 
tree  and  another  of  a  noncaterpillar  together  with  other  graphs  is  shown  in  Fig.  1.  It 
seems  that  Harary  and  Schwenk  were  among  the  first  to  study  these  trees  in  the 
mathematical  literature. ^'3,4 

In  chemistry  the  use  of  these  trees  resulted  from  studying  the  topological 
properties  of  benzenoid  hydrocarbons,  namely  resonance  relations  among  individual 
hexagons  of  a  benzenoid  system. 5  Two  hexagons,  in  a  benzenoid  hydrocarbon  are  called 
resonant  if  an  (aromatic)  sextet  (i.e.  a  set  of  three  circularly  conjugated  double  bonds) 
can  be  drawn  in  both  of  them  such  that  the  rest  of  the  carbon  atoms  are  spanned  either 
by  a  double  bond  or  by  a  sextet  of  electrons.  Gutman^  represented  such  resonance 
relations  among  hexagons  of  a  benzenoid  svstem  by  the  edges  of  a  caterpillar  tree: 
two  edges  in  a  caterpillar  tree  are  incident  if  and  only  if  the  corresponding  hexagons 
in  the  benzenoid  system  are  nonresonant.  Thus  the  tree  given  by  P4  (3, 0,4, 2)  corresponds 
to  the  benzenoid  hydrocarbon  B(P4(3,0,4,2))  drawn  in  Fig.  1.  There  is  a  one-to-one 
correspondence  between  the  labeling  of  the  edges  of  the  caterpillar  and  those  of  the 
hexagons  of  the  benzenoid  system.  Explicitly  these  terms  were  considered  in  chemistry 
(synonomously  under  the  name  "Gutman  trees")  in  three  recent  papers  by  this  author.^'® 

It  is  amazing  that  nearly  all  graphs  that  played  an  important  role  in  what  is  now 
called  "chemical  graph  theory"  may  be  related  to  caterpillar  trees.  For  this  reason 
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such  objects  are  of  great  importance  for  understanding  and  simplifying  combinatorial 
properties  of  much  more  complicated  graphs.  Three  main  areas  involve  uses  of  these 
trees,  viz.,  computational  methods,  ordering9  and  data  reduct  ion.  ^  it  is  convenient 
to  give  now  important  definitions  for  the  development  of  this  treatment. 


2.  Definitions  of  Important  Terms 
2.a  Graphs,  Lattices  and  Diagrams 


A  caterpillar  tree  might  be  associated  with  the  followi-g  objects: 


2.a.1  Polyhex  Graph  ^  B 


This  term  was  first  introduced  by  Hosoya  et  al.^  to  mean  a  graph  composed 
only  of  hexagons  which  have  an  even  number  2 1  of  points  and  thus  can  be  spanned 
by  disjoint  lines.  The  polyhex  B  is  one  such  type  while  B'  is  not  included  in  this 
definition 


Obviously  the  number  of  ways  in  which  B  can  be  spanned  by  the  £  disjoint  lines  is  well 
known  in  organic  chemistry  as  the  number  of  Kekul4  structures^  and  is  known  in 
mathematics  as  the  number  of  perfect  matchings,*^  a  synonym  for  a  famous  problem 
<n  dimer  statistics.^  Now  a  word  on  nomenclature: 

A  molecular  network  which  is  entirely  composed  of  hexagons  is  called  benzenoid. 
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If  no  three  hexagons  have  a  common  atom,  the  system  is  called  catacondensed.  If 
every  hexagon  of  a  catacondensed  system  has  at  most  two  neighboring  hexagons,  it 
is  said  to  be  nonbranched.  If  there  is  at  least  one  hexagon  in  a  catacondensed 
hydrocarbon  that  is  surrounded  by  three  other  hexagons  it  is  said  to  be  branched.  If 
in  a  polyhex  graph  at  least  one  vertex  is  common  to  three  hexagons  it  is  called 
pericondensed. 

2. a. 2  CJar_gra£hs^5,16 

Gutman^  seems  to  be  the  first  who  introduced  this  term  in  chemistry  and  the 
concept  was  later  developed  by  Gutman  and  this  author.^  For  nonbranched  benzenoids 
a  Clar  graph  is  simply  the  line  graph^of  a  caterpillar  tree.  In  fact  it  can  be  seen 
that  every  caterpillar  is  associated  with  a  Clar  graph  (c.f.  Fig.  1). 

For  branched  benzenoid  hydrocarbons,  however,  no  caterpillar  tree  is  defined, 
nevertheless,  a  Clar  graph  can  be  defined"^  in  the  following  way.  Let  h-|,  hj,  ...,  hn 
be  the  hexagons  of  the  branched  system.  Then  the  vertices  of  its  Clar  graph  are  v-|, 
v2«  •••»  vn  such  that  vj  is  connected  to  Vj  only  if  hj  and  hj  are  nonresonant. 

Relation  between  caterpillar  trees,  Clar  and  polyhex  graphs 

At  this  point  it  is  important  to  digress  on  the  relation  between  the  three  types 
of  graphs  defined  above.  We  start  by  considering  the  polyhex  graph  again.  First  we 
observe  that  two  hexagons  in  a  polyhex  may  or  may  not  be  resonant.  An  illustration 


is  considered  below 
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Thus  hexagons  3  and  5  at  left  are  nonresonant  while  1  and  4  at  right  are. 

Gutman  demonstrated^  the  following  fact:  Every  nonbranched  (catacondensed) 
polyhex  graph,  B,  which  contains  n  hexagons  is  associated  with  a  caterpillar  tree,  T, 
containing  n  +  1  vertices  such  that:  two  (or  more)  incident  edges  in  T  correspond  to 
two  (or  more)  nonresonant  hexagons  in  B  and  vice-versa,  i.e.  two  (or  more)  nonincident 
edges  in  T  correspond  to  two  (or  more)  resonant  hexagons  in  B.  For  example  consider 
B(P4(3,0,4,2))  of  Fig.  1  and  let  hj  refer  to  the  rth  hexagon.  One  observes  that  linearly 
fused  rings  cannot  be  resonant.  For  example  none  of  the  hexagons  in  each  of  the 
following  sets  can  be  resonant:  {  h-],  h2»  (13,  h4>;  {  114,  115);  {(15,  hg,  h 7,  hg,  hg,  h -]q } ; 
{h-jQ,  h-ji,  h-]2  )•  Similarly  we  say:  aJJ_  the  edges  in  T  in  each  of  the  following  sets 
are  adjacent:  { e-j,  e2»  eg,  04};  {e4,  eg};  {65,  eg,  ey,  eg,  eg,  e-jo);  el1<  e12^ 

This  one-to-one  correspondence  extends  to  the  vertices  of  the  Clar  graph:  All  the 
vertices  in  each  of  the  following  sets  are  adjacent  (refer  toA  (P4  (3,0, 4, 2))):  {  v-|,  V2» 
V3,  v4 };  {v4,  V5  };  {  V5,  vg,  V7,  vg,  vg,  v-^)}  )  (  v10'  v11’  v12  These  fundamental 
relations  have  important  implications  in  understanding  the  combinatorial  structures 
of  benzenoid  hydrocarbons  as  we  shall  see  later.  Because  of  this  relation  to  benzenoid 
systems,  caterpillar  trees  will  be  also  called  benzenoid  trees. 

2. a. 3  King  Polyomino  Graphs^®  p 

Consider  a  rectangular  lattice  composed  of  cells  arranged  in  certain  number 
of  rows  and  columns.  Such  graphs  are  called  polyominos  or  square  animals.  Two  cells 
in  a  polyomino  are  defined^®  to  be  adjacent  if  they  share  at  least  one  vertex.  The 
maximum  number  of  adjacent  cells  is  therefore  four.  This  corresponds  to  four 
nonresonant  hexagons  annellated  in  a  linear  fashion.  Because  of  this  fact  king 
polyominos  might  be  made  to  correspond  to  polyhex  graph  containing  linear  segments 
which  are  no  more  than  four-hexagons  long.  As  an  illustration  we  consider  the  following 
set  of  graphs  whose  caterpillar  tree  is  P4  (2,0,2, 1); 
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Therefore  for  every  subset  of  adjacent  cells  there  is  a  corresponding  subset  of 
nonresonant  hexagons,  incident  edges  and  adjacent  vertices. 


2.a.4  Rook  Boards^  pr 

A  rook  board  is  a  subset  of  cells  of  a  j  x  j  chessboard.  Godsil  and  Gutman 
demonstrated^  that  every  bipartite^  graph  G  is  associated  with  a  rook  board  such 
that  a  cell  Cjj  (which  is  located  in  the  ith  row  and  jth_  column  of  the  board)  exists  only 
if  vertices  i  and  j  are  connected  in  G.  We  illustrate  how  rook  boards  are  constructed 
which  correspond  to  bipartite  caterpillar  trees  in  the  following  chart 


1'  2  2  3 

T4(2, 1,0,1) 
V  3'  4 


*0*0 

1  2'  3  4' 

V  2'  2  4 


1'  2'  3'  4' 


In  the  above  Chart  three  different  labellings  of  the  bipartite  caterpillar  generated 
three  different  rook  boards.  These  are  shown  by  the  heavily  outlined  squares.  In  a 
rook  board  two  cells  are  adjacent  if  they  share  the  same  row  and  column.  For  example 
in  the  top  board  the  cell  labelled  3  is  adjacent  to  1,2, 4,5  while  cells  in  4  and  7  are 
not  adjacent.  The  three  rook  boards  preserve  the  combinatorial  counts  in  the  caterpillar 
form  which  they  were  generated.  Thus,  e.g.  there  are  six  triplets  of  nonadjacent  edges 
in  that  tree,  viz.,  {(146),  (147),  (157),  (257),  (247),  (246)}.  The  same  subsets  of  cells 
in  all  three  boards  are  nonadjacent.  Similarly  one  can  easily  demonstrate  that  there 
are  13  sets  of  nonadjacent  edges,  each  of  cardinality  2  in  the  tree  which  correspond 
to  13  such  sets  of  nonadjacent  cells  in  any  of  the  above  boards.  Such  combinatorial 
counts  of  nonadjacent  structures  represent  the  coefficients  of  counting  polynomials 
which  will  be  considered  later.  Idential  combinatorial  counts  exist  in  the  associated 
polyhex  graph,  Clar  graph  and  king  polyomino  graph,  all  shown  below 

2  4 

— o - o 

13  5  6  7 

B(T4  (2,1,0,!))  A(T4 (2, 1,0,1))  P  (T4 (2,1,0, !» 

There  are  one-to-one  correspondences  between  the  labellings  of  the  edges  of  T4  (2, 1,0,1), 
the  hexagons  of  B,  the  vertices  of  A  and  the  cells  of  Pr.  Thus  a  knowledge  of 
nonadjacent  edges  in  T  yields  details  of  nonadjacent  structures  in  other  graphs,  namely, 
polyhex  graph,  Clar  graph,  king  and  rook  polyominos. 

2.b  Polynomials  and  Nonadjacent  Structures 

All  polynomials  of  caterpillar  trees  and  related  graphs  and  lattices  are 
combinatorial  descriptors  of  the  "nonadjacent  structures"  in  a  given  object.  It  seems 
that  Hosoya^l  was  the  first  to  introduce  the  concept  of  a  nonadjacent  structure  in 
chemistry.  For  a  graph,  C,  he  defined  a  counting  polynomial,  H(G;x)  by 


,  4L  *»  I* 


.  »L,  .«  »'  '|L  «l.  «r  *»»  Alt'  tk’j't 


H(G;x)  =  ^  p(C;k)  xk 


where  p(G;k)  is  the  number  of  ways  of  selecting  k  nonadjacent  edges  in  G  (i.e.  k  edges 
in  which  no  two  of  them  are  adjacent).  The  term  p(G;k)  is  called  the  number  of 
k-matchings  in  G.  Conveniently  in  H(G;x)  (and  in  all  other  combinatorial  polynomials) 
p(G;o)  is  taken  to  be  unity  and  m  is  the  maximal  value  of  k.  A  more  general  expression 
of  such  polynomials^  is  given  by  Eq.  (2),  viz., 


max.  k 

F(G;x)  =  l  p  0(C;k)  x^k'n) 
k 


where  p,  0(G;k)  and  /(k,n)  are  all  functions  of  the  particular  polynomial.  Table  1  lists 
several  polynomials  of  some  use  in  chemistry.  As  an  illustration  we  consider  the  sextet 
polynomial^  Qf  a  benzenoid  system,  o(B;x).  This  important  polynomial  plays  quite 
a  significant  role  in  the  chemistry  of  benzenoid  hydrocarbons^  and  was  first  defined 
by  Hosoya  and  Yamaguchi^2  as  a  combinatorial  enumeration  of  the  number  of  Kekul4 
structures  of  a  benzenoid  hydrocarbon.  As  can  be  inferred  from  Table  1  and  Eq.  (2), 
d  =  1  and  /(k,n)  =  k  for  the  sextet  polynomial.  Further,  the  generating  function^4 
is  given  by  r(B,k)  xk  where  r(B,k)  is  called  the  k_th_  resonant  number  of  the  polyhex 
graph  of  the  benzenoid  system  B.  It  measures  the  number  of  selections  of  k  mutually 
resonant  and  disjoint  hexagons  in  B.  The  "nonadjacent  structures"  of  Table  1  may 
be  depicted  from  the  following  diagram 


Nonadjacent  Structures 


Edges 


Caterpillar 

(Benzenoid) 


Vertices 


Hexagons 


Benzenoid 


Cells 
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The  following  identities  can  easily  be  written 

H(T4(2,1,0,1);x)  *w  (A(T4(2,1,0,1);x))  = 
o  (B(T4(2,1,0,1);x))  =  K<Pr(T4(2/1,0,1);x))  = 

K(P(T4(2,1,0,1);x))  =  1  +  7x  +  13x2  +  6x3  . 

Naturally  when  we  set  x  =  1  in  the  generating  function  we  arrive  at  the  number  of 
Kekul£  structures,  K(B)  5  K,  a  problem  which  is  continuously  being  the  focus  of 
interest23  despite  its  early  history  in  chemical  combinatorics.2^  Thus  a  knowledge 
of  the  counting  polynomial  of  a  given  caterpillar  leads  to  other  polynomials  such  as 
sextet,  independence,  color,  king  and  rook  polynomials  (if  the  latter  two  boards  exist). 
The  above  treatment  which  applies  to  nonbranched  benzenoid  hydrocarbons  can  easily 
be  extended  to  other  systems  as  in  the  following. 

2.c  Branched  Systems 

By  the  application  of  the  appropriate  recursive  relations  of  the  sextet  polynomial23 
one  can  associate  a  "pseudobenzenoid"  tree  (i.e.  a  benzenoid  tree  containing  a  variable 
x)  with  virtually  any  benzenoid  hydrocarbon.  The  principle  is  simple:  choose  any  row 
of  hexagons  and  divide  the  set  of  Kekul6  patterns23  into  the  set  of  distinctive  cases 
so  that  each  vertical  line  in  that  row  is  chosen  double.  Caution  should  be  taken  for 
the  possibility  that  the  chosen  double  and  the  resultant  fixed  double  bonds  might  produce 
a  proper  sextet22  by  assigning  double  bonds  to  the  remaining  skeleton.  The  sextet 
polynomial  of  the  branched  benzenoid  hydrocarbon  can  be  written  in  terms  of 
polynomials  of  nonbranched  systems.  The  nonbranched  fragments  can  then  be 
transformed  into  caterpillar  trees  whose  counting  polynomials  are  identical  to  the 
sextet  polynomials  of  the  nonbranched  polyhex  graphs.  The  algorithm  is  illustrated 
in  Fig.  2  for  a  branched  system,  where  the  resulting  pseudobenzenoid  tree  is  shown 


in  braces. 


3.  Benzenoid  trees  and  Computation 


If  the  matching  polynomial  and/or  the  counting  polynomial  of  a  benzenoid  tree 
is  known  all  other  polynomials  of  related  graphs  in  Table  1  become  immediately 
available.  Now  the  matching  and  counting  polynomials  are  related  to  each  other  as^® 


a  (T;x)  =  xN  H(T,  -x‘“);  (3) 

KT;x)  =  i'N  xN/2  a(T;ix-j)  ,4, 


where  N  is  the  number  of  vertices  in  the  given  caterpillar  tree.  Using  elegant  operator 
algebra  defined  by  Hosoya  and  Ohkami23,28  lt  js  possible  to  compute  either  function*. 
(3)  and/or  (4)  for  families  of  caterpillar  trees.  As  an  illustration  we  calculate  tn*- 
matching  polynomials  of  the  family  of  caterpillars  of  Pn(3,3,  3)  (where  the  3's  a*-e 

repeated  n  times).  We  apply  the  recursion^  to  the  starred  edge  as  follows: 


(1) 


ttJ  =  •  •  VM...W  -  t  ■  m...w 
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Tn 
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n  n-l  1 
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Defining  the  step-up  operator  O  such  as 


O  ^n  =  ^n+1 


Steps  (1)-(4)  can  be  re-written  in  the  following  form: 


(O  +  x2)  Tn_-|  =  xKn  -  0  Jn  -  0  ln.-|  =  0 
x  Tn.-|  +  Kn  -  xJn  -  0  ln.-]  =  0 
Tn-i  +  0  Kn  +  Jn  -  x  ln.-|  =  0 
(x3  -  Ox)  Tn.-i  +  0Kn^0Jn  +  O  ln_i  *  0 


A  nontrivial  solution  of  (6)  requires  that 


A 

(O  + 


c2) 


-x 


0  0 


x  1  -x  0 

1  0  1  -x 

(x3-Ox)  0  0  0 


i.e. 


O2  ♦  (3x2  -  x4)0  +  x6  =  0 


Application  of  8  on  Tn  leads  to 


t.  ■  n  .*»  , 


'  *  f  ,  -•»  •’»  J  l.l  ‘  J  <4I  ‘.|  *.?  **»“(•  rn  i'a,  |*i  1*1.1  j.»  Ifc' I  >  i.1  i>*  »J’|J  |.|Vl\  t't  t'i 


Tn+2  ♦  <3x2  -  x4)Tn+1  ♦  x8Tn  =  0 


It  is  easy  to  show  that  for  the  general  case,  Tn  (m,m,  ...,  m),  the  operatar  eqn.  8, 
becomes: 


O2  +  (mxm_1  -  xm+1)  6  +  x2m  =  0 


Using  Eq.  (10)  and  the  following  two  identites: 


a(T -|(1);x)  =  x2  -  1;  a(T2(1,1);x)  =  x4  -  3x2  ♦  1 


then  repeated  application  of  eqn.  (10)  (for  m  =  1)  leads  to 


a(T8(1, 1,1,1, 1,1,1,1);x)  =  a(T8(l8>;  x)  = 


x16  .  15x14  +  74x12  .  290x10  +  258x8  -  290x6  ♦  74x4  -  15x2  ♦  1 


The  graphs  of  this  polynomial  are  shown  in  Fig.  3.  Using  relations  (3)  or  (4)  we  can 
write  the  following  identities  (See  Table  1): 

H(T8(18);x)  =u1(A(T8(18);x))  =  o(B(T8(18);x)) 

=  K(P(Tg(18);x))  =  K(Pr(T8(18);x)) 

=  1  +  15x  +74x2  +  290x3  +  258x4  +  290x5  +  74x8  +  15x7  +  x8  (12) 


So  we  know  immediately,  e.g.,  that  there  are  290  ways  of  placing  3  or  5  nonattack i 


«  *  \  »  _  «  _  a  _  •  .  •  .  •  .  •  -  »  .  »  A  * 


kings  on  the  chessboard  shown  in  Fig.  3.  It  is  interesting  to  observe  that  there  is  t 

only  one  way  of  selecting  a  subset  of  nonadjacent  structures  of  maximum  cardinality  | 

t 

(=  8  in  this  case.)  Polynomials  of  such  types  as  counting,  matching,  etc.  of  graphs  ; 

I 

are  symmetric  in  the  sense  of  Eq.  (13),  viz.,  ] 

t 

l 

I 

I 

aj  =  a^.j  (13) 

where  N  is  the  number  of  vertices  in  T  (or  N-1  =  the  number  of  hexagons  in  B  =  number 

I 

of  vertices  in  A  etc.).  Benzenoid  hydrocarbons  for  which  eqn.  (13)  holds  (i.e.  with  a 
"symmetric"  sextet  polynomial)  are  known  to  have  a  single  sextet  formula  (i.e.  a  single 
Clar  representation^).  For  such  types  the  number  of  aromatic  sextets  they  contain 
is  very  close  to  the  number  of  sextet-type  resonance  interactions  per  Kekul6  structure.^1 
This  last  statement  is  known  as  Aihara's  conjectureA^I  The  importance  of  Aihara's 
observation  is  because  it  specifies  a  condition  of  the  benzenoid  hydrocarbon  in  which 
case  the  simple  Clar  sextet  formalism-^  roughly  estimates  its  Dewar-type  resonance 
energy^  Gutman^  commented  on  Aihara's  conjecture  by  defining  a  function  F(B) 
by 


h 

F(B)  =  M(B)  K(B)  -  2  J  K(B  -  Hj)  (14) 

i 

where  M(B)  =  max.  k  =  the  maximum  cardinality  of  a  set  of  mutually  resonant  but 
disjoint  aromatic  sextets  in  the  benzenoid  graph  and  K(B)  is  its  Kekul4  count.  The 
summation  of  the  second  term  is  taken  over  all  hexagons,  Hj,  of  B  where  h  is  the  total 
number  of  Hj.  Gutman  specified  Aihara's  condition  that  the  hydrocarbon  be  represented 
by  a  single  sextet  formula  by  having 


r(B;  max.k)  =  1 


(15) 
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Fig.  4  shows  an  example  of  a  hydrocarbon  which  satisfies  Aihara's  conjecture  and  another 
that  does  not.  Furthermore,  Gutman  restated  Aihara's  postulate  by  defining  F(B)  = 
0  whenever  eqn.  (15)  holds.  Benzenoid  hydrocarbons  for  which  F(B)  vanishes  define 
the  "best  conditions"  where  nearly  all  existing  resonance-structure  theories  apply, 
viz.,  Clar's,30  Dewar's^2  conjugated  circuits  theory  of  Randid  (34)  as  well  as  the 
structure-resonance  theory  of  Herndon. ^5  In  fact  all  such  theories  coincide  for 
benzenoid  systems  for  which  Eq  (15)  applies. 

For  benzenoid  trees  which  do  not  possess  elements  of  symmetry  or  which  do 
not  belong  to  a  given  periodic  network  of  trees  the  method  of  Balasubramanian  and 
Randi6*6  becomes  particularly  suitable  for  computation  of  a(T;x).  For  example  for 
the  caterpillar  tree  P4  (2,0,2, 1)  shown  above  we  have  the  following  identities 


The  quotient  tree,  Q  =  P4 

“1  =  a3  =  ot(T i<2))  =  x^  -  2x  thus  a-|'  =  03'  =  x2 
a  2  =  (P'j(O))  =  x;  a2  =  1 

a  4  =  x2  -  1;  04'  =  x  (16) 


From  the  adjacency  matrix  of  Q,  a(T;x)  =  a  (P4<2,0,2,1);x)  is  given  by  the  following 
determinant 


-(x^  -  2x)  x2 

1  -x 

0  x2 

0  0 


0  0 

1  0 

-(x^  -  2x)  x2 

X  -(x2-1) 


(17) 


In  the  above  notation  the  primed  letter  denotes  the  matching  polynomial  of  a  type, 
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(such  as  1-1(2))  after  its  root  vertex  has  been  pruned.  The  notation  is  essentially  that 
used  in  ref.  36. 

4.  Caterpillar  Benzenoid  trees  and  the  Ordering  of  Graphs;  A  Relation  with  Young 

I 

Diagrams  | 

In  their  work  on  algebraic  characterization  of  skeletal  branching,  Gutman  and 
Randid^  used  theorems  of  Muirhead®®  to  order  and  compare  a  set  of  trees  (caterpillars 
and  noncaterpillars).  In  their  treatment  two  trees  are  characterized  by  a  sequence  | 

of  nonnegative  integers  {a-],  a2»  ...,  a^}  and  {b-|,  b2,  ...»  bj}  representing  the  degrees 
of  their  vertices  when  listed  in  descending  orders.  For  example  P4  (3, 0,4, 2)  would 
be  associated  with  the  sequence  {6,  4,  3,  2,  1,  1,  1,  1,  1,  1,  1,  1,  1}.  Muirheacfs  | 

conditions  state  that  Ta  whose  sequence  is  {a-|,  a2,  ...»  }is  greater  than  whose 
sequence  is  { b-j,  b2,  ...,  bj}  if 


a1  il  bi 

a1  +  a2  ^  bi  +  b2 


(18) 


ai  +  a2  +  ..  +  ak  =  t>1  +  +  ..  +  bj 

Whenever  am  +  an  +  ...  +  aQ  _>  bm  +  bn  +  ...  +  bQ  but  ar  +  as  +  ...  +  at  _>  br  +  bs  +  ... 
♦  bt,  the  two  tree  graphs  are  said  to  be  noncomparable.  The  latter  lead  to  bifurcation 
sites  in  the  ordering  hierarchy.  Using  the  above  criteria  Gutman  and  Randi6^  ordered 
sets  of  trees  for  which  N  =  8,9  and  10.  Furthermore  they®^  discovered  the  very 


interesting  observation  that  their  ordering  of  trees  can  be  made  to  overlap  with  Ruch 
and  Schdnhofer®^  ordering  of  a  set  of  Young  diagrams  if:  (a)  information  on  the  terminal 
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vertices  is  suppressed  and  (b)  the  valency  of  each  vertex  is  reduced  by  one.  This 
significant  result  leads  to  a  relation  between  a  Young  diagram  and  a  caterpillar  tree 
and  whence  between  a  Young  diagram  and  nearly  all  other  graphs  and  lattices  used 
in  chemistry  and  physics,  namely,  Clar  graphs,  king  polyomino  graphs,  rook  boards 
and  polyhex  graphs.  As  an  illustration,  the  Young  diagram  which  corresponds  to  the 
set  of  graphs  of  T  =  P4  (2,  0,  2,  1)  is  shown  below 


\LMS- 

(4,  3,  2,2,  1,1,1,  1,1,5 


(3,2,  I,  I.)  i 


which  will  be  denoted  as  Y  (3,2, 1,1).  There  is  a  unique  Young  diagram  for  every 
caterpillar  tree  (or  any  of  its  associated  graphs)  but  the  reverse  is  not  true,  i.e.  two 
(or  more)  caterpillars  may  be  related  to  the  same  Young  diagram.  The  following 
examples  (from  the  set  N  =  8)  illustrate  this: 

{  P3(4,0,1),  P3(1,3,1)}  £  Y(4,1,1)  ; 

{  P3  (3,  0,  2);  P3  (3,  1,  1),  P3  (2,  2,  1)}e  Y(3,  2,  1); 

{  P4(3,  0,  0,  1),  P4  (1,  2,  9,  1)}  e  Y(3,  1,  1,  1); 

{  P4(2,  0,  0  2),  P4(2,  0,  1,  1),  P4(2,  1,  0,  1),  P4(1,  1,  1,  1,)}  e  Y(2,  2,  1,  1); 

{  P5(2,  0,  0,  0,  1),  P5(1,  0,  1,  0,  1),  P5(1,  1,  0,  0,  1)}  c  Y(2,  1,  1,  1,  1) 

Fig.  5  shows  the  ordering  of  all  Young  diagrams  containing  six  boxes.  Fig.  6 
is  the  corresponding  order  of  the  nonbranched  benzenoid  systems  which  correspond 
to  the  caterpillar  trees.  The  numbers  in  parentheses  are,  respectively,  Y-|,  >2'  Y3 
and  >4,  where  Yj  is  a  permutation  integral  (of  Herndon-^)  involving  permutation  of 


(4i  +  2)  pi  electrons  in  the  benzenoid  systems.  Twice  these  numbers  enumerate  R*j, 
^2'  ^3  and  R4  (respectively)  where  Rj  is  a  conjugated  circuit^4  containing  (4i  +  2) 
pi  electrons.  It  is  obvious  that  the  second  number,  ^2  (or  R2)  is  almost  constant  for 
the  same  level  of  ordering  of  the  Young  diagram.  In  Fig.  6  the  polyhex  graphs  are 
represented  by  their  LA  sequences^.  Thus  a  hexagon  may  be  annellated  in  two  ways, 
viz., 


Linear  Angular 


and  by  convention^  the  terminal  hexagons  are  labelled  by  L.  Thus,  e.g.  B(P4<3,0,4,2)) 
is  denoted  by  LLLAALLLLALL  or  L^A^L^AL^.  In  Fig.  the  nonbranched  benzenoid 
hydrocarbons  containing  seven  hexagons  are  ordered.  The  numbers  in  parentheses 
are  respectively  >-],  Y  2,  Y3  and  T4  where  >j  is  a  Herndon  permutation  integral^ 
involving  permutation  of  (4i  +  2)  pi  electrons.  Naturally,  twice  these  numbers  lead 
to  the  corresponding  conugated  circuits, ^4  R-j,  R2,  R3  and  R4. 

It  is  emphasized  here  that  through  relating  benzenoid  trees  to  other  graphs  beside 
benzenoid  hydrocarbons,  such  as  Clar  graphs,  king  polyominos  and  rook  boards,  they 
can  all  be  ordered  according  to  schemes  of  Ruch  and  Schdnhofer  adopted  by  Young 

diagrams.-^ 

5.  Benzenoid  Trees  and  Data  Reduction 

An  important  part  of  an  analysis  of  chemical  data  is  the  data-reduction  step. 
In  the  past  this  involved  mainly  curve-fitting  procedures.  The  role  of  graph  theory 
was  recognized  in  the  work  of  Smolensk  ii4^  and  later  of  Cordon  and  Kennedy.4^  The 
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dualist  graphs  of  Balaban  and  Harary^2  might  also  be  regarded  as  a  type  of  structure 
reduction  representing  polyhex  graphs  of  benzenoid  hydrocarbons.  Recently  this 
author^  explored,  for  the  first  time,  the  possibility  of  using  benzenoid  trees  to  store 
and  retrieve  information  on  related  benzenoid  systems  (i.e.  a  bezenoid  system  whose 
sextet  polynomial  is  identical  to  the  counting  polynomial  of  the  tree).  Several  physical 
and  combinatorial  properties  including  electronic  absorption  spectra,  heats  of 
atomizations,  number  of  conjugated  circuits,  number  of  self-avoiding  walks,  number 
of  Sachs  graphs  are  studied  and  in  all  cases  excellent  correlations  are  found  between 
the  natural  logarithms  of  a  property  of  the  benzenoid  hydrocarbon  and  simple  powers 
of  the  connectivity  index  of  its  tree  graph. ^  As  an  illustration  Fig.  7  shows  a  plot 
of  the  number  of  Kekul4  structures  in  In  units  of  a  homologous  series  of  the  zigzag 
polyacenes  and  the  connectivity  indices  of  their  trees,  x(T)'s,  given  by 


x (T)  =  £(dj  djH 


where  the  summation  is  taken  over  all  edge  types  in  T,  (dj  dj)'s  (dj  is  the  degree,  i.e. 
valency  of  vertex  i  in  T). 


6.  Other  Avatars  of  Caterpillar  trees.  Generation  of  cycle  indices  of  wreath  product 

groups:^ 

The  composition^  of  two  groups  A  and  B  is  denoted  by  A [ B ]  (read:  A  around 
B)  is  known  as  the  "wreath  product"  or  the  "Cruppenkranz".  A  permutation  in  A[3] 
is  given  by 


<aj  ei,  62,  ...,  ScjXxj,  y,)  =  (axj,  6 j Y j) 


where  aeA,  Be  B  and  the  sequence  ...  may  not  involve  necessarily  distinct 
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elements.  The  elements  (xj;yj)  arise  from  the  cross  product  of  x  =  (  x-j,  X2,  ...,  xj) 
and  V  =  {  y),  y2»  ....  yeh  Wreath  product  groups  have  a  number  of  chemical  and  physical 
applications  discovered  recently  by  Balasubramanian^  who  revived  interest  in  Polya's 
Theorem. ^  We  show  here  that  caterpillar  trees  might  be  used  as  a  model  to  visualize 
the  cyclic  structure  and  the  operations  of  Eq.  (20).  we  take,  as  an  illustration  the 
group  S2  [S2 J.  Then  we  have:  x  =  (a,b) ,  Y  =  {1,2,} and 


x  x  Y  =  { al,  a2,  bl,  b2  }.  The  system  might  be  envisaged  as 


The  operations  in  S2's  are:  (1X2),  (12),  (a)(b),  (ab).  Then  the  following  2.2^  =  8  elements 


exist  in  S2tS2l,  viz.. 


(1)  ((a)(b);  (1X2),  (1X2)) 

(2)  «a)(b);  (12),  (12)) 

(3)  ((a)(b);  (1X2),  (12)) 

(4)  ((a)(b);  (12),  (1X2)) 


and  four  other  elements  using  the  operation  (ab)  instead  of  (aXb). 

The  element  (1)  is  simply  the  identiy  element  which  corresponds  to  the  operation 


Such  an  element  generates  six  one  cycles,  i.e,  s^.  The  element  (2)  operates  as  follows 
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((a)(b);  (12),  (12»  (a  1)  =  (a2) 
«a)(b);  (12),  (12))  (a2)  =  (al) 
((a)(b);  (12)  (12))  (bl)  =  (b2) 
((a)(b);  (12),  (12))  (b2)  =  (bl) 
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<(a)(b);  (12),  ('1X2))  (al)  =  (a2) 

((a)(b);  (12),  (1)(2))  (a2)  =  (al) 

((a)(b);  (12),  (1)(2))(b1)  =  (bl) 

((a)(b);  (12),  (1)(2))(b2)  =  (b2) 

Hence  this  element  is  represented  by  (al  a2)(b1)(b2)  and  also  contributes  s-|2  S2-  The 
caterpillar  model  is  shown  below 


a  b  a  b 


This  "caterpillar-modeling"  of  the  operations  of  the  composition  of  two  groups  facilitates 
the  understanding  of  the  abstract  algebra  involved  in  the  definition  especially  for 
beginners.  A  similar  model  is  considered  by  Balasubramanian^  where  he  uses  a 
"particle-in-box"  model,  thus: 


12  2  1 


a  b  A  B 


Actually  either  model  generates  the  permutation  group  of  the  non-rigid  N2H4  molecule^ 
(i.e.,  the  nitrogen  atoms  are  represented  by  the  root  of  vertices  of  P2<2,2)  while  its 
monovalent  vertices  represent  the  hydrogen  atoms.).  The  above  modeling  can  be 
extended  as  shown  below 


mm  m  m  m 
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Extension  to  the  recently  defined^  generalized  wreath  product  is  also  possible.  Thus 
the  NMR  group  of  butane  might  be  modeled  by  P4(3,2,2,3)  which  represents  S2  [S3, 

S21. 

7.  Conclusion 


Although  caterpillar  (Benzenoid  =  Gutman)  trees  are  not  widely  known  in  the 
chemical  literature  its  uses  span  a  wide  range  of  applications  including  data  reduction, 
computations,  ordering  and  modeling  notations  of  abstract  groups  (such  as  wreath 
and  generalized  wreath  product  groups)  which  are  necessary  for  NMR  spectrscopy 
and  counting  distereomers.49 
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Table  1  Some  polynomials,  associated  graphs  and  graph  invariants  in  chemistry  and  physics 


Table  1  (cont.) 


a  Number  of  selections  of  k  independent  edges  e  T  (i.e.  no  two  edges  are  incident) 

b  Number  of  selections  of  k  nonadjacent  but  mutally  resonant  hexagens  B 

c  Number  of  ways  of  arranging  k  nontaking  kings  on  a  polyomino  graph 

d  Number  of  ways  of  arranging  k  non-attacking  kings 

e  Number  of  selections  of  k  independent  vertices  eA  (No  two  are  adjacent) 
f  Number  of  colorings  in  C  in  which  there  are  k  vertices  of  the  same  color  so 
that  no  two  of  them  are  adjacent. 
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A  caterpillar  tree,  P 4(3,0,4,2),  the  corresponding  Clar  graph,  A  (P4O,  0,  4,  2))  and 
the  corresponding  benzenoid  hydrocarbon,  B(P4(3,0,4,2)).  A  noncaterpillar  tree  is  also 
shown. 

Fig.  2 

Recursive  generation  of  "pseudocaterpillar"  tree  of  a  branched  benzenoid  hydrocarbon. 
The  factor  of  X  accounts  of  the  proper  sextet^^  jn  the  graph  to  the  right. 

Fig.  3 

Illustration  of  eqn.  (13)  for  a  caterpillar  tree  (i.e.  benzenoid  tree)  and  its  associated 
graphs  for  Max  k  =  8.  Simple  application  of  eqn.  (10)  shows  that  there  are  290  ways 
of  placing  either  3  or  5  non-attacking  kings  on  the  chessboard  P  or  Pr.  The  subset 
of  invariants  leading  to  X8  is  heavily  outlined. 

Fig.  4 

Examples  of  benzenoid  hydrocarbons  which  possess  one  Clar  representation  (1  and 
2)  and  a  hydrocarbon  with  two  Clar  representations  (3a  and  3b).  Hydrocarbons  1  and 
2  satisfies  Aihara's  conjecture  (eqn.  15). 

Fig.  5 

Ruch's  ordering  of  all  Young  diagrams  containing  six  boxes.  Site  of  bifurcations  indicate 
noncomparable  diagrams. 

Fig.  6 

Ordering  of  nonbranched  benzenoid  hydrocarbons  which  are  in  one-to-one  correspondence 
with  the  Young  diagrams  shown  in  Fig.  5.  The  polyhex  graphs  are  denoted  by  their 
L-A  sequences.8  Numbers  in  parentheses  are  (Yi,  Y 2>  Y4>  respectively,  where 

i  is  a  Herndon  permutation  intregral88  involving  permutation  of  (4i  +  2)  pi  electrons. 


Twice  these  numbers  lead  to  (R-],  R2,  R3,  R4):  the  sequences  of  the  corresponding 


conjugation  circuits. 


34 


Fig.  7 

Correlation  between  In  K(B)  i.e.  the  natural  logarithms  of  the  Kekul4  counts  of 
zigzag  polyacenes  (1  =  phenanthrene,  2  =  chrysene,  3  =  picene,  4  =  fulminene,  ... 
X(T):  the  connectivity  indices  of  the  relevant  caterpillar  (i.e.  benzenoid  trees). 
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